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A 40 cm ion thruster is being developed at the NASA Glenn Research Center to obtain 
input power and propellant throughput capabilities of 10 kW and 550 kg, respectively. The 
technical approach here is a continuation of the “derating” technique used for the NSTAR 
ion thruster. The 40 cm ion thruster presently utilizes the NSTAR ion optics aperture 
geometry to take advantage of the large database of lifetime and performance data already 
available. Dome-shaped grids were chosen for the design of the 40 cm ion optics because 
this design is naturally suited for large-area ion optics. Ion extraction capabilities and 
electron backstreaming limits for the 40 cm ion optics were estimated by utilizing NSTAR 
30 cm ion optics data. A preliminary service life assessment showed that the propellant 
throughput goal of 550 kg of xenon may be possible with molybdenum 40 cm ion optics. 
One 40 cm ion optics’ set has been successfully fabricated to date. Additional ion optics’ sets 
are presently being fabricated. Preliminary performance tests were conducted on a 
laboratory model 40 cm ion thruster. 


Nomenclature 

Adome Domed grid area, m 2 

d Grid aperture diameter, m 

d h Beam extraction diameter, m 

e Electron charge, C 

Eeff Effective elastic modulus of the perforate 

grid region. Pa 

f corT Sputter yield correction factor 

Fdome Fractional change of the grid area from 
doming 

ffl at Beam flatness parameter (ratio of average- 

to-peak beam current density) 

Accelerator grid view factor 
Grid open area fraction 
Beam divergence thrust correction factor 
Acceleration of gravity, m/s" 

Height of the grid dome, m 
J b Beam current, A 

j b Beam current density. A/m" 

Ji Accelerator grid impingement current, A 

k B Boltzmann's constant, J/K 

kcorr Correction factor for peak beam current 

density, V- m 2 /A 

l e Effective acceleration length, m 


l g Grid cold gap, m 

l n Neutralization length, m 

l up Upstream potential well location, m 

niamu Atomic mass of grid material, amu 
m d Discharge chamber mass flow rate, kg/s 
mi Ion mass, kg 

m n Propellant atomic mass, kg 

N Neutral density distribution function 

n Number of gravities (i.e., g-load) 

N a Avogadro’s number, mol” 1 

n b Facility background neutral density, m“ 3 

n d Discharge chamber neutral density, nf 1 

n n Neutral density distribution upstream of the 

neutralization plane, irr 

NPH Normalized perveance per hole, A/V " 
NPHcorr Correlated normalized perveance per hole, 
A/V 3/2 

r c Chord radius of domed grid, m 

Rdome Radius of curvature for the grid dome, m 

Rj Ratio of doubly-to-singly charged ion 

current 

R°max Maximum net-to-total accelerating voltage 

extrapolated to a zero beam current 
T Thrust, N 

T g Grid temperature, K 
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T n Neutral propellant temperature, K 

T s Stiffener temperature, K 

t s Screen grid thickness, m 

V b Beam voltage, V 

Va Discharge voltage, V 

V eb Electron backstreaming limit voltage, V 

V t Iim Impingement-limited total voltage, V 

V'_, Accelerator grid pit and groove volumetric 
erosion rate per aperture, mVs 
V' Screen grid volumetric erosion rate per 

aperture, mVs 

V', Accelerator grid aperture wall volumetric 
erosion rate per aperture, mVs 
Y Sputter yield at an orthogonal angle of 

incidence 

z Axial dimension, m 


a Thrust correction factor for doubly-to-singly 

charged ion current 

a g Grid material thermal expansion coefficient, 

1/C° 

a s Stiffener material thermal expansion 

coefficient, 1/C° 

rid** Discharge chamber propellant utilization 

efficiency accounting for doubly charged 
ions 

©max Maximum dome angle from the grid axial 
centerline, rad 

v e ff Effective Poisson's ratio of the perforate 
grid region 

n Pi (3.141593) 

p Grid material mass density, kg/m 3 

peff Effective mass density of the perforate grid 
region, kg/m 3 

a c . e Charge-exchange cross-section, m“ 

ahoop” Peak circumferential (i.e., hoop) stress, Pa 
CTmend^ Peak meridional stress. Pa 
<p Neutral transparency of the ion optics 

(pi Screen grid transparency to discharge ions 


Superscripts 

Accelerator grid 
avg Average 

d n Charge-exchange ions created downstream 

of the neutralization plane 
h >! Temperature during thruster operation 


Orthogonal to the grid surface 
NSTAR ion thruster 


Peak 


n Room temperature 

s Screen grid 

up Charge-exchange ions created upstream of 

the neutralization plane 
Per aperture 

+ Singly charged ion 

Doubly charged ion 

Subscripts 

a Accelerator grid 

s Screen grid 

Introduction 

The success of the NSTAR (i.e., NASA Solar Electric 
Propulsion Technology Applications Readiness 
Program) 30 cm ion thruster system on the Deep 
Space 1 mission has demonstrated the viability of ion 
propulsion for deep space missions. 1,2 As a result, ion 
propulsion is a candidate for several deep space 
missions, such as the Neptune Orbiter, Titan Explorer, 
Mars Sample Return, Europa Lander, and others. 
However, ion propulsion system mass and volume 
could be significantly reduced for many of these 
missions by increasing the NSTAR thruster's 
propellant throughput and peak input power 
capabilities beyond the 140 kg and 2.3 kW, 
respectively, demonstrated by past and ongoing wear 
tests. 3-5 Propulsion system mass and volume 
reductions occur because fewer thrusters and, 
therefore, fewer accompanying power processors and 
propellant feed system components, would be 
required. 

Increasing propellant throughput and thruster power is 
limited, in part, by sputter erosion of the ion optics. 6,7 
Significant charge-exchange accelerator grid sputter 
erosion can lead to electron backstreaming due to 
accelerator aperture enlargement, grid structural 
failure due to pit and groove erosion of the 
downstream surface, or an unclearable grid short by a 
flake from sputter-eroded accelerator grid material. 
Significant screen grid sputter erosion by discharge 
chamber ions can lead to structural failure of the 
screen grid, as well as exacerbating the 
aforementioned accelerator grid failure mechanisms 
from direct beam ion impingement. 

Long duration tests in the NSTAR program 3-5 and 
others 6,7 have demonstrated, however, that sputter 
erosion of the downstream surface and aperture walls 
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of the accelerator grid must be addressed to increase 
ion optics' longevity. Wear test results from Ref. [3] 
have shown that screen grid sputter erosion is minimal 
for the NSTAR thruster, likely due to the derated 
operating conditions of this engine. Although screen 
grid erosion must be addressed in any ion optics' 
service life assessment, lifetimes well beyond those 
demonstrated in Refs. [3] to [5] are expected as long as 
the low discharge voltages and plasma densities are 
maintained and the ratio of doubly-to-singly charged 
ions is not significantly increased. 

Several methods of extending ion optics' service life 
are presently being pursued at the NASA Glenn 
Research Center (GRC). 8 One method involves 
changing the ion optics' material to one with a lower 
volumetric sputter erosion rate than molybdenum. Grid 
material technologies presently investigated include 
carbon-based materials, titanium, and ion-implanted 
materials. 8 " 11 Another method for increasing ion optics 
service life involves increasing the accelerator grid 
thickness to provide more material for sputter 
erosion. 1- This prevents the loss of the large database 
already available with molybdenum ion optics while 
providing an anticipated 2x improvement in 
accelerator grid service life. While all of these 
methods allow for increases in thruster propellant 
throughput and input power, most improvements are 
anticipated to be modest if the 30 cm thruster diameter 
is maintained. 9 

To obtain input power and propellant throughput 
capabilities of 10 kW and 550 kg, respectively, a 
40 cm ion thruster is being developed at NASA 
Glenn. 13 This engine has twice the beam area as the 
NSTAR 30 cm ion thruster, whose beam diameter is 
28 cm. The technical approach here is a continuation 
of the “derating” technique used for the NSTAR ion 
thruster. With regards to the ion optics, this approach 
maintains low beam current densities across the ion 
optics to extend ion optics service life. Therefore, the 
40 cm ion thruster can, at minimum, double input 
power and propellant throughput per thruster for the 
same thrust densities as the NSTAR thruster. The 
40 cm thruster can further offer increased input power 
and propellant throughput when other grid 
technologies are used, such as thick-accelerator-grid 
(TAG) ion optics. 12 

This paper reports on the preliminary results of the 
development of 40 cm ion optics for the 40 cm ion 


thruster. Results to date include ion optics design 
considerations, and preliminary service life 
assessments and test results. 

Design Considerations 

Ion Thruster Performance Goals 

The design goals of the 40 cm ion thruster include a 
1 to 10 kW input power throttling range and a 550 kg 
xenon propellant throughput capability. Emphasis was 
placed on providing a flat beam profile to further 
improve ion optics service life. 1 ' This improvement 
occurs because several ion optics failure mechanisms 
are a result of the local sputter erosion of the ion optics 
that are a function of the local beam current density. 

Ion Optics’ Geometry Performance 

The 40 cm ion thruster is presently envisioned to use 
the NSTAR ion optics aperture geometry, described in 
Refs. [14] and [15], and the thick-accelerator-grid 
geometry, described in Ref. [12]. The following 
analyses will focus on the use of the NSTAR 
geometry. Analyses with this aperture geometry allow' 
the 40 cm thruster to take advantage of the large 
database of lifetime and performance data already 
available with this geometry. 

The ion extraction capability of the 40 cm ion optics 
was analyzed first. Because the ion optics' aperture 
geometry of the NSTAR thruster was utilized, 30 cm 
ion optics data could be used to predict 40 cm ion 
optics' performance. Average beam currents per 
aperture (or beamlet currents) as a function of the 
impingement-limited total voltage for an NSTAR ion 
thruster are shown in Fig. I. 111- The impingement- 
limited total voltages were defined as the voltage 
where the slope of the accelerator current as a function 
of total voltage (i.e., the sum of the absolute values of 
the beam and accelerator power supply voltages) was 
-0.02 mA/V. 


The ion extraction capability of a given aperture 
geometry can be correlated to other operating 
conditions by utilizing a normalized perveance per 
hole (or NPH): 


NPH = 



3 


ylim -* 



( 1 ) 


where is the beamlet current, is the 

impingement-limited total voltage, d s is the screen 
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aperture diameter, and l e is the effective acceleration 
length. The effective acceleration length is defined as: 



where 1« is the grid gap and t* is the screen grid 
thickness. The impingement-limited normalized 
perveance per hole for the NSTAR ion optics is plotted 
in Fig. 2 as a function of discharge-to-total voltage 
ratio. As the figure demonstrates, the normalized 
perveance per hole is a function of the discharge-to- 
total voltage ratio. This effect has been noted in other 
studies 16,17 and can be compensated by defining a 
correlated normalized perveance per hole (or NPHc 0rT ): 


NPH 


corr 


= NPH 


( 


\ 


A. 

wlim 

V t 


V* 

j 


(3) 


where V d is the discharge voltage. 16 In Ref. [16], q was 
set to 0.2, however, a value of 1.2 was found to yield 
constant correlated value in this study. The correlated 
normalized perveance per hole for the NSTAR ion 
optics is also plotted in Fig. 2. As Fig. 2 shows, the 
correlated normalized perveance per hole is 
approximately 2.4x1 0“ n A/V 3/2 for the NSTAR ion 
optics" aperture geometry. 


Total voltages for the 40 cm thruster may be 
determined by using the correlated normalized 
perveance per hole to determine the impingement- 
limited total voltage for a given beam current, and then 
adding voltage margin to the predicted limit. One 
selection criterion for the total voltage at a given beam 
current was that the total voltages provide a minimum 
perveance margin (i.e., the difference between the total 
voltage during normal operation and the impingement- 
limited total voltage) of at least 200 V. 


High specific impulses are typically best suited for the 
deep space missions anticipated for the 40 cm ion 
thruster. This generally results in thruster operation at 
high R-ratios (i.e., net-to-total accelerating voltages), 
and, therefore, high accelerator voltages. These high 
accelerator voltages are beneficial because they reduce 
accelerator grid erosion rates and, therefore, extend grid 
service life. The upper limit of the accelerator voltage is 
limited, however, by the onset of electron 
backstreaming. This electron backstreaming limit can 
be estimated with a semi-empirical equation derived by 
Kaufman for the maximum R-ratio as a function of the 
ion optics" aperture geometry. 18 However, since the 


40 cm ion optics' aperture geometry is the same as that 
for the NSTAR ion optics, electron backstreaming limit 
data from NSTAR ion optics' tests can be used. 11,12 


The electron backstreaming limit and, therefore, the 
maximum R-ratio, are also functions of the peak 
beamlet current. This peak beamlet current is directly 
proportional to the peak beam current density just 
downstream of the ion optics, and, therefore, 
proportional to the beam current density divided by the 
beam flatness parameter (i.e., average-to-peak beam 
current density). 10,11 Data with the NSTAR 30 cm ion 
optics from Ref. [11] shows that this relation is a 
function of beam power supply voltage and may be 
nonlinear. However, the electron backstreaming can be 
conservatively estimated by assuming a -1 V per 
1 A/m“ dependence on the peak beam current density. 
The electron backstreaming limit (or V e b) for the 
40 cm ion thruster can, therefore, be estimated by: 


v e „ = v b . 


( \ D 0 
1 ” K rmx 


R° 


4-J, 


n ' fflat 


■•k„ 


(4) 


Here, Jb is the beam current, db is the beam extraction 
diameter, fn at is the beam flatness parameter, k^ OIT is the 
correction for peak beam current density (i.e., -1 V per 
1 A/nr), and R 0 ^ is a maximum R-ratio extrapolated 
to a zero beam current. The value for R 0 ^ for the 
molybdenum 30 cm ion optics tested in Ref. [11] was 
determined to be 0.918. Equation (4) predicts 30 cm 
ion optics' electron backstreaming limits to within 
+ 15V/-5V. 


For operation at high R-ratios, the selected accelerator 
voltage was determined by adding voltage margin to 
the predicted electron backstreaming limit. 

A preliminary throttle table for the 40 cm ion thruster 
was developed based, in part, on the predicted 
performance of the ion optics' aperture geometry. 13 
Results relevant to the ion optics are shown in Table 1 
for operation at high R-ratios. Assumptions for these 
proposed throttling points include a 24 V discharge 
voltage and a beam flatness parameter of 0.6. As the 
table shows, the lowest perveance margin was 
estimated to be 209 V. This value occurred at the 
lowest thruster input power. The lowest electron 
backstreaming margin (i.e., the difference between the 
absolute values of the accelerator voltage and electron 
backstreaming limit) occurred at the highest beam 
voltage and current and was estimated to be 40 V. 
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Grid Design Considerations 

Dome-shaped grids, similar to those used on the 
NSTAR thruster, were chosen for the design of the 
40 cm ion optics. This was because dome-shaped grids 
are naturally suited for large-area ion optics. Reasons 
for this included: 1) the dome shape mitigates 

thermally induced buckling during thruster operation, 
and, therefore, provides a relatively stable grid gap 
throughout the ion optics' area; and 2) the dome shape 
provides stiffness under vibration loads during launch. 
Furthermore, there is design heritage with this grid 
design at NASA Glenn. 


As discussed in previous sections, the NSTAR ion 
optics' aperture geometry was utilized for the 40 cm 
ion optics. As a result, the beam extraction diameter, 
grid thicknesses, and perforated geometries were fixed 
for the 40 cm grids. The only remaining, adjustable 
parameter was the dome height. It will be shown in the 
following paragraphs that selecting a dome height-to- 
chord radius similar to that of the NSTAR grids 
addressed most key issues associated with grid 
fabrication and thruster operation while having little 
impact on other issues concerning thermal expansion 
and launch stresses. 


Fabrication Issues 

The grid dome shape is hydrostatically formed, a 
process developed by NASA Glenn in 1972, 19 and the 
apertures are photo-chemically etched. The primary 
concern regarding grid fabrication was that over- 
forming the 40 cm grids could result in aperture 
misalignment, significant residual stresses, and 
possibly even fracture the grid material. Since a grid is 
under a biaxial state of stress during forming (grid 
thicknesses are small enough to assume a thin-walled 
membrane), then the fractional change in area as a 
result of the forming process is the relevant parameter 
here. It can be shown that the area of the domed grid 
(or Aaome) is given by: 

A dome =rc-(h 2 +r;), (5) 


where h is the dome height and r c is the radius at the 
dome base, or the chord radius. The fractional change 
in area from forming (or F d ome) is, therefore, given by: 


1 dome 


71 (h 2 +r 2 )-7t*r 2 


7t ■ r : 



(6) 


As the equation shows, the fractional change in 
area from forming is a function of the dome 


height-to-chord radius. Maintaining the same ratio as 
for the NSTAR 30 cm grids will, therefore, mitigate 
over-forming the 40 cm grids. 


Thruster Operation Issues 

Two thruster operation issues were identified for the 
40 cm ion optics. These included thrust losses due to 
the dome shape and thermally-induced changes in the 
grid gap. 


The thrust losses considered here are caused by the 
shape of the dome, which directs beamlets off-axis. 
The following analysis is similar to that done in 
Ref. [19]. The apertures of the NSTAR and 40 cm 
screen and accelerator grids are designed to be aligned 
throughout the grids' radius, with no beamlet steering 
to collimate the beam. As a result, beamlets are 
directed normal to the surface of the grids. The thrust 
(or T) produced by a domed grid can, therefore, be 
expressed as: 


T = 

(7) 


Here, V b is the beam voltage, mi is the ion mass and e 
is the electron charge. Within the double integral, 0^ 
is the maximum angle from the centerline of the grids, 
j b n is the beam current density normal to the grid 
surface, Rdome is the radius of curvature of the dome, 
a is the double-to-single ion current correction factor, 
and F t n is the beamlet divergence correction factor that 
does not account for divergence from the grid dome 
shape. The latter term is used to separate beam 
divergence due to beamlet expansion from beam 
divergence due to the dome shape. To simplify this 
analysis, it will be assumed that the beam current 
density normal to the grid surface and the ratio of 
double-to-single ion current are constant throughout 
the grid. The beam current density can, therefore, be 
reduced to the beam current divided by the area of the 
dome. Equation (7) can be expressed as: 

T = EKiK f; ft- Rj ome -sin 2 e mx . (8) 

V e Adome 


It can be shown that the last two terms in the equation 
above reduce to r c \ Using this and equation (5), 
equation (8) reduces to: 


J 2 ' v ‘ 

1 + 

r h ' 

2 

If e 





( 9 ) 
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As this equation shows, the thrust is a function of the 
dome height-to-chord radius. Maintaining the same 
ratio as that for the NSTAR 30 cm grids will, 
therefore, keep these thrust losses due to the dome 
shape of the grids the same for the 40 cm grids. 


Grids thermally expand during thruster operation, 
which can change the grid gap when the grids are 
heated to their final operating temperatures. An 
analysis was conducted to compare 30 cm and 40 cm 
grids under thermal loads. The following simplifying 
assumptions were made (similar to those in Ref. [19]): 
1 ) each grid and its mounting ring were at uniform, but 
different, temperatures; 2) the final shape of the heated 
grids is spherical; 3) both grid diameters investigated 
had the same temperatures and temperature 
distributions; and 4) there are no thermally-induced 
meridional stresses. The chord radius and arc length 
across the grids (or S) during operation are given by: 

i;'” = r" • (l + <x s • AT S ), and (10) 


S hot =S" (l + a g AT g ). (11) 

Here, a s and a g are the thermal expansion coefficients 
for the stiffener and grid materials, respectively, while 
AT S and AT g are the change in temperatures for the 
stiffeners and grids, respectively. The “it” and “hot" 
superscripts denote parameter values at room 
temperature and temperatures during thruster 
operation, respectively. With equations (10) and (11), 
it can be shown that: 


f r c n - (l + ct, • AT„ ) ) 

11 D hrt I K dome 

k •'dome ) 


= arc sin 




( 12 ) 

where the radius of curvature for the heated grid is 
solved iteratively. Since it is known that: 

h — ^ dome _ ^ dome ~ " *c" » (13) 


the height of the heated grid can be solved with 
equations (12) and (13). 


Temperatures for this analysis were taken from 
Ref. [20] for a NASA J -Series 30 cm ion optics set 
operated at a discharge power of 210 W, which is 
similar to the NSTAR thruster at a 1.44 kW input 
power of and a 1.1 A beam current. Radial grid 
temperatures profiles were curve-fit with second order 
polynomial equations and integrated to determine 
average grid temperature changes. The thermal 
expansion of the grid-mounting hardware was also 


included in the analysis. Note that the results for the 
40 cm grids would not be applicable to the same 
discharge power as for the 30 cm thruster, but would 
correspond to a much higher discharge power (and, 
therefore, beam current) since the 40 cm grids have a 
larger thermal mass. Furthermore, the temperature 
distribution across the grids and mounting hardware 
would not necessarily be the same because of the 
different thermal environments of the two thruster 
sizes. Regardless, the aforementioned assumptions are 
only intended to provide a comparison for the two grid 
diameters. 


Analysis results are shown in Fig. 3, where the fraction 
of the hot-to-room temperature grid gap is plotted as a 
function of dome height-to-chord radius. The predicted 
fraction for the 30 cm ion optics was at the lower 
uncertainty limit of the measured value, indicating that 
the model reasonably predicts thermally-induced grid 
gap changes at this power level. As Fig. 3 further 
demonstrates, the predicted results for the 40 cm ion 
optics show that: 1) there is very little change in the 
fraction of the hot-to-room temperature grid gap over a 
wide range of dome height-to-chord radius ratios; and 
2) the overall hot grid gap is about 2% smaller than 
that for the 30 cm ion optics at steady-state conditions. 
Because the dome height-to-chord radius ratio has 
little effect on the grid gap during steady state thruster 
operation, it is reasonable to utilize the same ratio as 
that for the NSTAR 30 cm grids so that fabrication 
issues and thrust losses are mitigated. 


Launch Stress Issues 

The two spacecraft launch issues examined included 
vibration-induced stresses and grid movement. In 
order to simplify the analysis, it was assumed that final 
shape of the grids under these stresses was spherical, 
which requires grid deflections to be small. 
Furthermore, only grid vibrations in the axial direction 
were examined. Since the apertures significantly 
weaken the perforated grid regions, these regions were 
analyzed. The perforated regions were assumed to 
have a tangential top edge support. For thin-walled 
spherical domes under axial loads from there own 
mass, peak meridional and circumferential (i.e., hoop) 
stresses (or Omend*™* and Choop™ 1 *, respectively) are 
given by: 21 


max 
me rid 


rc ’ § P eff ^ dome 

i + cos(e niax ) 


, and 


(14) 
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max _ A ’ g * P e ff * ^ dome 
^ hoop — ^ 


(15) 


Here, n is the number of gravities (i.e., g-load), g is the 
acceleration of gravity, and p e n is the effective mass 
density of the grid perforated region. Note that 
l+cosCBmax) ~2 for all cases examined here, so that it 
can be assumed that the meridional and 
circumferential stresses are approximately equal. 
Solving equation (13) for R^ome, It can be shown that: 


1 + 


'h v 


max max 

mend ~~ u hoop 


; • P eff 


• r„ 


(16) 


r c 


As the above equation shows, maintaining a constant 
dome height-to-chord radius will still result in a stress 
increase proportional to the chord radius. An 
unchanged dome height-to-chord radius would 
increase meridional and circumferential stresses by 
about 40% for the 40 cm grids. However, increasing 
dome height-to-chord radius to account for the 
increased stresses in the 40 cm grids would worsen 
thruster performance. To determine the significance of 
these increased stresses, maximum meridional and 
circumferential stresses were calculated for these grids 
as a function of g-loading. These values were then 
used to determine the maximum webbing stresses" 
and compared to the yield strength of the grid material 
to determine if the 40 cm grids were at risk of plastic 
deformation under launch loads in the axial direction. 
The resulting maximum stresses were found to be 
significantly smaller than the material yield stress, 
even at high g-loads. 


Note that this analysis omitted stresses resulting from 
radial motion, which should be included in any 
complete analysis. Unfortunately, this was beyond the 
scope of this initial effort. 


Vibration-induced axial motion was also examined. 
The same assumptions as the prior analysis were 
made. For thin-walled spherical domes under axial 
loads from their own weight, the change in dome 
height (or Ah) is given by: 21 


Ah = 


n g P eff -RL 


s«n 2 (0 max )+(l + v eff ) ln 


l + cos^x) j 


(17) 


Here, E e n and v e ff are the elastic modulus and 
Poisson's ratio, respectively, for the perforated region 
of the grids." 2 It can be shown that: 


1- 


'h' 


cos 


( 0 =- 


1 + 




v r c ) 


— . and 


S * n 2 ( 0 m ax)= 1 - 


1- 


'h V 


1 + 




(18) 


(19) 


Both equations are merely functions of the initial 
dome height-to-chord radius. Solving equation (13) for 
Rdome. equation (17) can be given: 


Ah = 


n *g*P 


eff 


■ * r.. 


"eff 



( 20 ) 


where: 


k = sin : (0 nm )+(l + v eff )• In 


1 + cos(0 max ) 


.( 21 ) 


As equations (18) to (21) show, maintaining a constant 
dome height-to-chord radius will still result in a 
change in dome height proportional to the chord radius 
squared. Maintaining a constant dome height-to-chord 
radius would, therefore, approximately double grid 
movement for the 40 cm grids. If this grid movement 
is large enough, the grids could make contact under 
vibration loads. This occurrence, however, is not 
necessarily detrimental to the grids. Both the NSTAR 
30 cm grids and 50 cm were found to make contact 
under vibration loads without detrimental effects. 23,24 
As a result, this is not expected to be an issue for the 
40 cm grids. 


Preliminary Service Life Assessment 


A preliminary service life assessment was conducted 
to determine the lifetime gain by utilizing 40 cm ion 
optics. Service life analyses included sputter erosion of 
the downstream surface (i.e., pits and grooves) and 
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aperture walls of the accelerator grid, shown in Fig. 4, 
and sputter erosion of the upstream surface of the 
screen grid. 


Accelerator Grid 

Accelerator grid erosion is caused by charge-exchange 
ions created both upstream and downstream of the 
neutralization plane that bombard the accelerator grid. 
Accelerator grid erosion from direct impingement of 
beam ions is assumed to be negligible here. Charge- 
exchange ions created upstream of the neutralization 
plane sputter-erode both the aperture walls and the 
downstream surface of the accelerator grid, while 
those created downstream of the neutralization plane 
only sputter-erode the downstream surface (see Fig. 4). 
This assertion has been demonstrated, in part, by 
several numerical models.' 5 "" 8 It was also verified via 
testing where a two-grid system that was tested at high 
background pressures exhibited excessive pit and 
groove erosion of the downstream accelerator grid but 
little accelerator aperture enlargement. 29,50 


The volumetric erosion rate in the region of a single 
accelerator grid aperture (or V' ccel ) can be expressed 
as: 


V' = (v' up + v' up )+ v' dwn 

* accel V vv T p-g / T T p-g * 


( 22 ) 


Here, and V'_ g are the volumetric erosion rates of 

the aperture walls and pit and groove erosion, 
respectively, in the region of an aperture. The “up” and 
“dwn” superscripts denote that charge-exchange ions 
were created upstream and downstream of the 
neutralization plane, respectively. The variables in this 
equation can be expressed as: 


v /up + V /up = 

v w ^ v p-g 
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w/dwn 

V P"g 


y'dwn 
* i 

e-pa 


m amu y x* p-g 

n a ' ’ corr 


(24) 


other sputtering phenomena that can affect volumetric 
sputter erosion. Note that the prime superscripts 
denote parameter values per aperture. The 
aforementioned equations could be expressed without 
the prime superscripts to define the volumetric erosion 
rates for the entire accelerator grid. 


All of the above variables are known except for the 
impingement currents per aperture and the sputter 
yield correction factors. These impingement currents 
are from charge -exchange ions created both upstream 
and downstream of the neutralization plane. These 
currents will be determined here as was done by 
Monheiser and Wilbur, 5152 but with some 
modifications. The charge-exchange ion production 
current per aperture is given by: 

In 

J' up = J b a c _ e • jn n dz , and (25) 

Kp 

j; dwn =j;<wjF g n b dz. (26) 

Here, i\ is the beamlet current, a c _e is the charge- 
exchange cross-section, n n is the neutral density 
distribution upstream of the neutralization plane from 
discharge chamber neutrals, n h is the vacuum facility 
background neutral density, z is the axial dimension, 
lup is the upstream location of the potential well that 
traps charge-exchange ions in this region, and l n is the 
neutralization length. F g is the accelerator grid view 
factor for a charge-exchange ion created downstream 
of the neutralization, and, therefore, represents the 
probability that a charge-exchange created in this 
region will impinge the accelerator grid. The 
assumptions for equation (26) include, in part: 

1 ) potential gradients in the downstream beam plasma 
are too small to divert charge-exchange ions; and 

2) the facility background neutral density has isotropic 
velocity distribution. 


Here, J'i is the impingement current per aperture, 
is the atomic mass of the accelerator grid material, p a 
is the mass density of the accelerator grid, e is the 
electronic charge, N A is Avogadro’s number, Y is the 
grid sputter yield for an orthogonal angle of incidence, 
and f corT is a sputter yield correction factor. This 
correction factor accounts for non-orthogonal angle- 
of-incidence sputtering, variations in ion sputtering 
energies, re-deposition of sputtered material, and any 


Equations (25) and (26) differ from those derived by 
Monheiser and Wilbur in the following ways. First, the 
integral in equation (25) has l up as the upstream limit 
and not the accelerator grid downstream surface. This 
is because this upstream axial location more accurately 
defines those charge-exchange ions trapped in this 
region by the electrode potential gradients. Second, the 
integral in equation (25) does not include the 
contribution of the facility background neutral density. 
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This simplifying assumption is valid if background 
facility pressures are kept low. This appeared to be the 
case for the NSTAR thruster operated at NASA 
Glenn's Vacuum Facility 5, where comparisons of 
accelerator currents produced in this facility to those 
produced on the Deep Space 1 spacecraft were within 
0.4 mA (e.g., less than 10% of the total impingement 

at 2 kW). 2 Finally, the contribution of neutrals 

escaping the discharge chamber downstream of the 
neutralization plane were omitted in equation (26). 
This is because these neutrals do not have an isotropic 
velocity distribution, but a velocity distribution 

directed away from the accelerator grids, and, 

therefore, should be omitted. 


It will be assumed here, as was done by Monheiser and 
Wilbur/ 1,32 that the neutral density distribution along 
the axis of an aperture can be approximated by the 
discharge chamber neutral density, n d , multiplied by a 
distribution function, N(z), so that: 

i n i n 

jn n -dz = n d • J N(z)-dz , (27) 

kp kp 

where N(z) is only a function of the ion optics' 
aperture geometry. Furthermore, the same one- 
dimensional view factor utilized by Monheiser and 
Wilbur will be used here. Equations (25) and (26) can, 
therefore, be re-written as: 

i n 

j; up = j; • O c _e • n d ■ | N(z) dz , and (28) 


j /dwn __ j b ~ ^c-e ’ ^ b 
1 ~ 2 

where d b is the beam extraction diameter. It can be 
shown that l n « d b , so that equation (29) can be 
expressed as: 

j/dwn __ J b ^c-e * n b ’ (30) 

1 4 

Because equation (30) utilizes a one-dimensional view 
factor approximation, it will only accurately predict 
impingement current contributions for apertures at the 
accelerator grid center, assuming small beamlet 
divergence half-angles. Equation (30) over-predicts 
impingement current contributions for apertures 
radially away from the accelerator grid center. 
Equation (30) also over-predicts impingement currents 
because potential gradients in the beam plasma were 



assumed to be too small to divert charge-exchange 
ions from the accelerator grid. Plasma potential 
measurements in NSTAR thrusters have demonstrated 
that there are radial and axial potential gradients/ 33 
Although these voltage gradients are only on the order 
of a few volts, charge-exchange ions produced in this 
region have significantly lower energies and can, 
therefore, be directed away from the grids. 

The discharge chamber neutral density is determined 
by: 


(>- 

-T|d + )m d 

<p-<k 

Jt m n • k B ■ T n 

32 


Here, is the discharge chamber propellant 

utilization efficiency accounting for doubly charged 
ions, m d is the discharge chamber flow rate, <p is the 
ion optics' neutral transparency, mn is propellant 
atomic mass, k B is Boltzmann's constant, and T n is the 
neutral propellant temperature. 


Because charge-exchange production is a direct 
function of beamlet current, the peak beamlet current 
is necessary for determining the maximum local 
sputter erosion rate of the accelerator grid. This peak 
beamlet current is given by: 



(32) 


where d a is the accelerator aperture diameter, J b is the 
beam current, and F oa a is the open area fraction of the 
accelerator grid. 


Equations (31) and (32) can be combined with 
equation (28) to yield: 


j'up-pk _ _a_ 


(l-nr)n 


f • F a 

1 flat 1 02 






m n k n • T n 


J N(z) dz ’ 


32 


(33) 

This equation can be used with equation (23) to define 
this peak volumetric sputter erosion rate per aperture. 
This is because the ion flux from this equation will 
impinge the aperture region within which it was 
created. However, the impingement current in 
equation (30) can impinge on any pit-and-groove site 
on the accelerator grid. It is, therefore, more 
convenient to utilize an average beamlet current for 
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equation (30). The resulting impingement current can 
be expressed as: 


j/dwn-avg _ * ^c-c 


* b *^a 


4 * d K • F“ 


(34) 


The effect of beam extraction diameter on accelerator 
grid volumetric erosion rates can now be assessed. 
Comparisons are made at the same operating electrode 
voltages, beam currents, and discharge chamber flow 
rates, so that the relative effects of beam extraction 
diameter can be determined. Because the ion optics' 
aperture geometry is unchanged, a peak volumetric 
erosion rate per aperture ratio can be established for 
charge-exchange ions produced upstream of the 
neutralization plane, where: 


v' up + v' up 

w T p-g 


(v- + v-r“ 


/ KICT*D 



c™ (i-ir) 

' (j_ nr - T -)- 


(35) 

Here, the “NSTAR” superscripts denote parameters for 
the NSTAR 30 cm ion thruster. Note that 
equation (35) assumes that changes in 1^ and l n do not 
affect the neutral density distribution integral in 
equation (28). It can be shown that this is likely the 
case for l n due to the long neutralization plane lengths 
and the exponential decay of neutral density 
downstream of an accelerator aperture. Although 
comparisons are made for the same ion optics' 
geometries and electrode voltages, it is presently 
unclear what the exact effect the varying beamlet 
current (and, therefore, varying space-charge) will 
have on l up . 


An average volumetric erosion rate per aperture ratio 
can further be established for charge-exchange ions 
produced downstream of the neutralization plane, 
where: 

w/dwn , NSTAR 

= (36) 

., dwn NSTAR j ' 

V P-S 

It is assumed that background neutral densities are 
independent of beam extraction diameter, which is 
reasonable because discharge chamber flow rates are 
assumed unchanged. 


The results of equations (35) and (36) are shown in 
Fig. 5. In this figure, the discharge chamber propellant 
utilization efficiencies and flatness parameters were 
assumed to be the same for all beam extraction 


diameters. As Fig. 5 demonstrates, increasing the 
NSTAR thruster beam extraction diameter to 40 cm 
reduces peak volumetric erosion rates by charge- 
exchange ions originating upstream of the 
neutralization plane by a factor of 0.25, or by 75%. 
This is expected since the number of accelerator 
apertures available for erosion was doubled and the 
discharge chamber neutral density available for 
charge-exchange production was halved. 

The 40 cm beam diameter further reduces average 
volumetric erosion rates by charge-exchange ions 
originating downstream of the neutralization plane by 
a factor of 0.7, or about 30%. This occurs because 
although the downstream accelerator grid surface area 
available for sputter erosion doubled for the 40 cm 
thruster, its view factor for downstream charge- 
exchange ions increased by a factor directly 
proportional to the 40 cm-to-30 cm ion optics diameter 
ratio. 

The service life enhancement due these reduced 
erosion rates requires knowledge of which 
impingement current (i.e., equation (28) or (30)) is 
more dominant. It may be possible to gain insight into 
the relative magnitude of each impingement current 
through short-term testing, which is the topic of 
Appendix A. However, a rough estimate of the 
improvement in accelerator grid propellant throughput 
capability may be made for the 40 cm ion optics if it is 
assumed that operating conditions are the same as 
those of the NSTAR thruster. Two accelerator grid 
failure mechanisms are considered separately. 

The first failure mechanism considered is electron 
backstreaming due to accelerator grid aperture 
enlargement. Since an ongoing life test of an NSTAR 
ion thruster has successfully processed about 140 kg of 
xenon, it will be conservatively assumed that this 
140 kg is the propellant throughput capability of the 
NSTAR ion optics with regards to electron 
backstreaming. Upstream of the neutralization plane, it 
is assumed that the fraction of charge-exchange ions 
that sputter-erode the aperture walls (as opposed to 
those that erode the downstream surface) is either 
constant or decreases with a decreasing beamlet 
current. This assumption is reasonable since NSTAR 
wear test results have shown that aperture wall mass 
loss due to sputter erosion exhibited a smaller flatness 
parameter (i.e., average mass loss per unit grid area 
divided by the peak mass loss per unit grid area) than 
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the beam flatness parameter. 3 ' 14 As a result, the ratio of 
accelerator aperture wall volumetric erosion rate per 
hole of the 40 cm thruster to that of the NSTAR 
thruster is equal to equation (35). Assuming that 
discharge chamber propellant utilization efficiencies 
and beam flatness parameters are the same for both 
thrusters, the 40 cm thruster offers a 4x increase in 
propellant throughput. The propellant throughput 
capability of the 40 cm thruster is, therefore, estimated 
to be about 560 kg for this failure mechanism, which 
exceeds the goal of 550 kg. 

The second failure mechanism considered is structural 
failure of the grid due to the pit and groove erosion of 
the downstream surface of the accelerator grid. The 
failure criterion here is assumed to be groove erosion 
through 80% of the grid thickness. Pit erosion is not 
considered because erosion through the grid thickness 
does not compromise the structural integrity of the 
grid. Past analyses based on NSTAR thruster wear test 
results estimated a propellant throughput capability of 
the NSTAR thruster of about 410 kg for this failure 
mechanism. 9 If it is conservatively assumed that 
charge-exchange ions created downstream of the 
accelerator grid dominate groove erosion, then the 
40 cm offers a 1 .4x increase in propellant throughput 
capability. The propellant throughput capability of the 
40 cm thruster is, therefore, estimated to be about 
570 kg for this failure mechanism, which exceeds the 
goal of 550 kg. 

The propellant throughput capability of the 40 cm 
thruster is estimated to be about 560 kg for electron 
backstreaming and 570 kg for a structural failure. The 
40 cm thruster is, therefore, anticipated to meet the 
propellant throughput goal of 550 kg of xenon. (Note 
that these preliminary results are rough estimates 
subject to the assumptions made for the calculations.) 
A more accurate service life assessment requires that 
detailed analyses be made at the anticipated operating 
conditions for the 40 cm ion thruster. 

Screen Grid 

Screen grid erosion is caused by discharge chamber 
ions that sputter erode the upstream surface of the 
screen grid. This sputter erosion generally causes a 
chamfering of the screen grid apertures. Failure occurs 
when the screen grid erosion results in the removal of 
the screen grid webbing, allowing beam ions to 
directly impinge on the accelerator grid. NSTAR 


thruster life testing indicated that screen grid erosion 
was peaked at the center on the grid/ This is expected 
since the center of the grid is where the radial beam 
current density, and discharge ion current density just 
upstream of the ion optics, is highest. 


Discharge chamber particles that sputter-erode the 
screen grid include both singly and doubly charged 
ions. The screen grid is at discharge cathode potential 
and the discharge chamber plasma is typically 
assumed to be at the discharge anode potential. As a 
result, singly charged ions are assumed to have 
energies equivalent to the discharge voltage while 
doubly charged ions have twice the discharge voltage. 
The volumetric erosion rate in the region of a single 
screen grid aperture (or V' ) can be defined as: 


v; = - 


e-ps 


m an 

N , 


• y t: 


s 

, "i, 

2 e p s ' N ; 


fS++ 

corr 


(37) 


Here, J s is the screen grid current per aperture, nT*™ 
is the atomic mass of the screen grid material, p s is the 
mass density of the screen grid, and f/orr is a sputter 
yield correction factor, similar to that used in the 
previous section. The “+” and “++" superscripts 
denote parameter values for the singly and doubly 
charged discharge chamber ions, respectively. 


Note that NASA xenon ion thrusters are typically 
operated at discharge voltages of about 25 V, which is 
a part of the derated approach. As a result, singly 
charged ions have energies on the order of the sputter 
yield threshold energies. 35 Furthermore, sputter yields 
increase exponentially with ion energy at these 
low discharge voltages. Because the ratio of doubly- 
to-singly charged ion current ranges between 0.05 to 
0.25 for the NSTAR ring-cusp discharge chamber,' it 
is reasonable to assume that doubly charged ions 
dominate screen grid erosion for these thrusters. 
Equation (37) can, therefore, be rewritten as: 


V>. 


rn 


2 e p s N, 


(38) 


Because screen grid erosion is highest at the center of 
the grid, a peak current is necessary for determining 
the maximum local sputter erosion rate of the screen 
grid in equation (38). It can be shown that this peak, 
doubly charged screen grid current per aperture is 
given by: 
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(39) 


Here, (ft is the screen grid transparency to discharge 
ions, d, is the screen grid aperture diameter, F s oa is the 
screen grid open area fraction, and Rj is the ratio of 
doubly-to-singly charged ion current. Note that 
equation (39) assumes that Rj is constant throughout 
the beam. Equation (39) can be plugged into 
equation (38) to determine a peak screen grid 
volumetric erosion rate per aperture. 


The effect of beam extraction diameter on the peak 
screen grid volumetric erosion rate can now be 
assessed. Comparisons are made at the same operating 
discharge voltages, total voltages, and beam currents 
so that the relative effects of beam extraction diameter 
can be determined. Because the ion optics aperture 
geometry is unchanged, a peak screen grid volumetric 
erosion rate per aperture ratio can be established: 
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The results of equation (40) are shown in Fig. 6. In this 
figure, the screen grid discharge ion transparencies, 
ratio of doubly-to-singly charged ion currents, and 
flatness parameters were assumed to be constant for all 
beam extraction diameters. The screen grid discharge 
ion transparency assumption is conservative because 
this ion transparency is known to increase with 
decreasing beam current densities, which can be 
inferred from data in Refs. [11] and [12]. As Fig. 6 
demonstrates, increasing the NSTAR thruster beam 
extraction diameter to 40 cm reduces peak screen grid 
volumetric erosion rates by a factor of 0.5, or by 50%. 
This is expected since the amount of screen grid 
webbing available for erosion was doubled. 

The service life of the 40 cm thruster is assumed to be 
double that of the NSTAR screen grid, given the 
aforementioned assumptions. NSTAR wear test results 
showed that screen grid erosion was minimal after 
processing 88 kg of xenon. 3 From these results, the 
NSTAR screen grid was estimated to have a 460 kg 
propellant throughput capability utilizing a 


conservative failure criterion of sputter erosion 
through half of the screen grid thickness. 9 The 40 cm 
thruster is, therefore, anticipated to be capable of 
processing about 920 kg of xenon, which exceeds the 
goal of 550 kg. 


Preliminary Results 

Fabrication of 40 cm Ion Optics 

One 40 cm ion optics’ set has been successfully 
fabricated to date, and is shown in Fig. 7 with 30 cm 
ion optics for comparison. The dome height-to-chord 
radius was matched to that of the NSTAR 30 cm ion 
optics. Additional ion optics’ sets are presently being 
fabricated. These additional optics’ sets not only 
include the NSTAR aperture geometry, but also a 
thick-accelerator-grid (i.e., TAG) aperture geometry 
for further improvements in accelerator grid service 
life. 12 

Preliminary Tests 

Test Hardware 

Preliminary performance tests were conducted on a 
laboratory model 40 cm ion thruster, described in 
Ref. [13]. Because 40 cm ion optics were unavailable 
for these initial tests, 50 cm ion optics from a prior 
NASA program were used. 24 These ion optics had an 
ion optics’ aperture geometry similar to the NSTAR 
geometry, with two differences. First, the grid gap was 
15% larger than the nominal NSTAR gap at the grid 
center and 8 to 33% larger at mid-radius. Second, the 
accelerator grid thickness was 25% thinner than the 
nominal NSTAR grid thickness. 

A power console similar to that described in Ref. [36] 
powered the thruster. This power console was 
modified to allow the thruster to be throttled up to 
10 kW. During the tests reported herein, neutralizer 
common was electrically tied to facility ground. A 
high purity gas feed system was used to provide xenon 
to the discharge cathode, discharge chamber, and 
neutralizer through separate mass flow controllers. 
Testing was conducted in Vacuum Facility 11 at 
NASA Glenn. This 2.2 m diameter by 7.9 m long 
facility is evacuated with seven cryogenic pumps and a 
turbomolecular pump. The total measured facility 
pumping speed was greater than 100,000 1/s with 
xenon. 


N AS A/TM — 200 1-211275 


12 



During thruster operation, voltages and currents were 
measured with digital multimeters and xenon flows 
with mass flow meters. Beam current density profiles 
were measured with a planar probe mounted onto a 
two-axis probe motion system. Both the probe and 
two-axis probe motion system are described in detail 
in Refs. [11] and [12]. The planar probe of this study 
had an exposed surface area of 0.25 cm". 

Test Results 

The two measurements reported below include 
impingement-limited total voltages and radial beam 
current density profiles. Although accelerator currents 
and electron backstreaming limits were measured, high 
background pressures, the larger-than-nominal cold 
grid gap, and the smaller-than-nominal accelerator grid 
thickness made comparisons with NSTAR ion optics’ 
data difficult. Impingement-limited total voltages were 
determined from plots of accelerator current as a 
function of total voltage where the slope was 
-0.02 mA/V. Unceitainties in impingement-limited 
total voltage determinations (and, therefore, perveance 
margins) were within ±10 V. All radial beam current 
density profiles were measured 49 mm downstream of 
the geometric center of the accelerator grid. 

Impingement-limited total voltages were measured at 
beam currents of 1.20, 2.70, 3.52, and 4.75 A. 
Although impingement-limited total voltages for this 
ion optics' set were decreasing with accumulated 
operation, there was no attempt to bum-in this optics’ 
set prior to collecting the data presented in this 
paper. 11,1 " As a result, perveance data are conservative 
and, therefore, preliminary. The resulting average 
beamlet currents as a function of the impingement- 
limited total voltage are plotted in Fig. 8 with NSTAR 
ion optics' data for comparison. As Fig. 8 shows, 
impingement-limited total voltages for the 40 cm 
thruster were 150 to 210 V higher than those of the 
NSTAR thruster ion optics. This is expected given the 
larger cold grid gap of the ion optics used on the 
40 cm thruster. Although the accelerator grid used on 
the 40 cm thruster was thinner than that of the NSTAR 
thruster, accelerator grid thickness has been shown to 
have a negligible effect on perveance. 12 ’ 37 38 

To permit comparisons between these 40 cm thruster 
data and the NSTAR ion optics' data, the 40 cm 
thruster impingement-limited total voltages were 
corrected for an NSTAR grid gap. This was done by 
using equation ( 1 ) to determine a normalized 


perveance per hole at a given beam current, and then 
using this value to determine a new impingement- 
limited total voltage for an effective acceleration 
length equal to that of the NSTAR ion optics. These 
corrected results are also plotted in Fig. 8. As the 
figure shows, corrected impingement-limited total 
voltages for the 40 cm thruster were less than 50 V 
higher than those of the NSTAR thruster ion optics. 

The mean value of the correlated normalized 
perveance per hole for the ion optics of the 40 cm 
thruster was 1.8x1 O' 11 A/V v2 , which is 25% lower 
than the 2.4xl0 -11 A/V v2 measured for the 30 cm 
NSTAR ion optics. The cause for this difference is 
presently unknown. Regardless, the 40 cm thruster 
with the 50 cm ion optics was able to successfully 
operate at the lowest total voltages in Table 1 for all 
beam currents tested, albeit with a reduced perveance 
margin. 

Radial beam current density profiles at beam currents 
of 1.20, 2.70, 3.52, and 4.75 A are shown in Fig. 9. 
Beam power supply and accelerator voltages were 
those corresponding to the minimum total voltages for 
each beam current listed in Table 1. Although the 
beam current density profiles exhibited a repeatable, 
slightly irregular profile within a 150 mm radius, there 
was no substantial peak in the profile center as with 
the NSTAR 30 cm ion thruster. 10-12 

Conclusions 

A 40 cm ion thruster is being developed at NASA 
Glenn to obtain input power and propellant throughput 
capabilities of 10 kW and 550 kg, respectively. This 
engine has twice the beam area as the NSTAR 
thruster, whose beam diameter is 28 cm. The technical 
approach here is a continuation of the “derating” 
technique used for the NSTAR ion thruster. With 
regards to the ion optics, this approach maintains low 
beam current densities across the ion optics to extend 
ion optics service life. The 40 cm thruster can further 
offer increased input power and propellant throughput 
when other grid technologies are used, such as thick- 
accelerator-grid (TAG) ion optics. 

The 40 cm ion thruster utilizes the NSTAR ion optics 
aperture geometry to take advantage of the large 
database of lifetime and performance data already 
available. The ion extraction capability and electron 
backstreaming limits of the 40 cm ion optics were, 
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therefore, be estimated by utilizing NSTAR 30 cm ion 
optics data. 

A preliminary throttle table for the 40 cm ion thruster 
was developed based, in part, on the predicted 
performance of the ion optics" aperture geometry. The 
lowest perveance margin was estimated to be 209 V 
and occurred at the lowest thruster input power. The 
lowest electron backstreaming margin was estimated 
to be 40 V and occurred at the highest beam voltage 
and current. 

Dome-shaped grids were chosen for the design of the 
40 cm ion optics because they are naturally suited for 
large-area ion optics. A dome height-to-chord radius 
ratio similar to that of the NSTAR grids was further 
utilized because it addressed most key issues 
associated with grid fabrication and thruster operation 
while having little impact on other issues concerning 
thermal expansion and launch stresses. 

A preliminary service life assessment was conducted 
to determine the lifetime gain by utilizing larger area 
ion optics. Utilizing long duration test results with 


NSTAR ion optics, comparisons indicated that the 
propellant throughput goal of 550 kg of xenon is 
possible with 40 cm ion optics utilizing the NSTAR 
ion optics geometry. 

One 40 cm ion optics" set has been successfully 
fabricated to date. Additional ion optics" sets are 
presently being fabricated. 

Preliminary performance tests were conducted on a 
laboratory model 40 cm ion thruster. Because 40 cm 
ion optics were not available for these initial tests, 
50 cm ion optics were used. These ion optics had a 
cold gap that was 8 to 33% larger than that of the 
NSTAR ion optics geometry. Impingement-limited 
total voltages were measured at beam currents of 1.20, 
2.70, 3.52, and 4.75 A. The 40 cm thruster with the 
50 cm ion optics was able to successfully operate at 
the lowest predicted total voltages for all beam 
currents tested. Impingement-limited total voltages 
corrected for their larger cold grid gap for the 40 cm 
thruster were less than 50 V higher than those of the 
NSTAR thruster ion optics. 
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Appendix A 


Accurately assessing the accelerator grid service life 
enhancement with 40 cm ion optics requires knowing 
whether charge-exchange ions created upstream or 
downstream (i.e., equation (28) or (30), respectively) 
of the neutralization plane dominate the total 
impingement current. It may be possible to gain 
insight into the relative magnitude of each 
impingement current by comparing impingement 
currents of both 30 and 40 cm ion optics at similar 
operating conditions. Total impingement currents can 
be determined for equations (28) and (30) by replacing 
the beamlet current with the beam current, to yield: 

T U P i ~ (l “*ld )* 

J i J b °c-e j = 
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Because the N(z) and l up are unknown, solving for the 
value of each current is not possible. However, 
impingement current ratios can be used to determine 


the relative magnitude of each term as a function of 
beam extraction diameter. Using the same assumptions 
as with equations (35) and (36), it can be shown that: 
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Equations (A3) and (A4) show that if 40 cm thruster 
accelerator currents are greater than those for the 
NSTAR thruster at similar operating conditions 
(assuming similar flatness parameters and discharge 
propellant utilization efficiencies), then charge- 
exchange ions originating downstream of the 
neutralization plane likely dominate the total 
impingement current. However, if 40 cm thruster 
accelerator currents are approximately half those for 
the NSTAR thruster at similar operating conditions, 
then charge-exchange ions originating upstream of the 
neutralization plane dominate the total impingement 
current. 
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Table 1. Preliminary 40 cm thruster ion optics’ currents and voltages 
for operation at high R-ratios.“ 


Input 
Power, b 
kW 

Beam 

Current, b 

A 

Beam 

Voltage, b 

V 

Accelerator 
Voltage, 11 V 

Total 

Voltage, 11 

V 

Total 
V oltage 
Limit, c V 

Electron 
Backstreaming 
Limit, d V 

10.0 

5.80 

1567 

-257 

1824 

1196 

-217 

8.2 

4.75 

1567 

-257 

1824 

1111 

-203 

6.1 

3.52 

1567 

- 257 

1824 

994 

-186 

4.7 

2.70 

1567 

-257 

1824 

901 

-175 

2,2 

1.20 

1567 

- 257 

1824 

667 

-156 

7.4 

4.75 

1396 

-243 

1639 

1111 

-187 

5.5 

3.52 

1396 

-243 

1639 

994 

-171 

4.2 

2.70 

1396 

-243 

1639 

901 

-160 

2.0 

1.20 

1396 

-243 

1639 

667 

-140 

4.7 

3.52 

1179 

-227 

1406 

994 

-152 

3.6 

2.70 

1 179 

-227 

1406 

901 

-141 

1.7 

1.20 

1179 

-227 

1406 

667 

-121 

3.2 

2.70 

1021 

-216 

1237 

901 

-127 

1.5 

1.20 

1021 

-216 

1237 

667 

-107 

1.1 

1.20 

679 

-197 

876 

667 

-76 


a A discharge voltage of 24 V and a flatness parameter of 0.6 were assumed. 
^Preliminary operating point. 

Predicted impingement-limited total voltage based on the average beamlet current. 
Predicted limit. 
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Fig. 1. Average beamlet current as a function of impingement-limited total voltage 
for the NSTAR ion optics. Data taken from Refs. [11] and [12]. 
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Fig. 2. Normalized perveance per hole as a function of 
discharge-to-total voltage ratio for the NSTAR ion optics. 
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Fig. 3. Fraction of hot-to-room temperature grid gap for 30 cm and 40 cm ion optics. 

Measured datum is from Ref. [20]. 
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Ratio of Accelerator Grid Volumetric 
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Fig. 6. Ratio of peak screen grid volumetric erosion rate per hole to that of the NSTAR thruster. 
Screen grid discharge ion transparencies, ratios of doubly-to-singly charged ion currents, discharge 
voltages, and flatness parameters are assumed to be the same for all beam extraction diameters. 



Fig. 7. Photograph of 40 cm screen grid (left foreground) and 40 cm accelerator grid (right foreground) 
with a 30 cm accelerator grid (center background) for comparison. 
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Fig, 8. Average beamlet current as a function of impingement-limited total voltage 
for the 40 cm and NSTAR ion thrusters. Data for the 40 cm thruster include 
measured data and data corrected for an NSTAR effective acceleration length. 



Radius, mm 


Fig. 9. Radial beam current density profiles for the 40 cm ion thruster at 49 mm downstream 
of the ion optics’ center for beam currents of 1.20 A, 2.70 A, 3.52 A, and 4.75 A. 
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